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ABSTRACT 

The interaction of a rotating star’s magnetic field with a surrounding plasma disk lies at the heart of 
many questions posed by neutron stars in X-ray binaries. We consider the opening of stellar magnetic 
flux due to differential rotation along field lines coupling the star and disk, using a simple model for the 
disk-opened flux, the torques exerted on the star by the magnetosphere, and the power extracted by 
the electromagnetic wind. We examine the conditions under which the system enters an equilibrium 
spin state, in which the accretion torque is instantaneously balanced by the pulsar wind torque alone. 
For magnetic moments, spin frequencies, and accretion rates relevant to accreting millisecond pulsars, 
the spin-down torque from this enhanced pulsar wind can be substantially larger than that predicted 
by existing models of the disk-magnetosphere interaction, and is in principle capable of maintaining 
spin equilibrium at frequencies less than 1 kHz. We speculate that this mechanism may account for 
the non-detection of frequency increases during outbursts of SAX J1808.4-3658 and XTE J18I4-338, 
and may be generally responsible for preventing spin-up to sub-millisecond periods. If the pulsar wind 
is collimated by the surrounding environment, the resulting jet can satisfy the power requirements 
of the highly relativistic outflows from Cir X-1 and Sco X-I. In this framework, the jet power scales 
relatively weakly with accretion rate, Lj oc and would be suppressed at high accretion rates 

only if the stellar magnetic moment is sufficiently low. 

Subject headings: magnetic fields — accretion, accretion disks — pulsars: general — stars: neutron 
— stars: jets — X-rays: binaries 


1. INTRODUCTION 

The study of accreting pulsars was inaugurated by the 
discovery of 4.8 s X-ray pulsations from Cen X-3 ( Giac- 
coni et al. 1971). These first observations immediately 
revealed sigmhcant pulse period changes on timescales 
as short as hours and a secular decrease in period, char¬ 
acteristics that would become central to understanding 
this new source class. Shortl y thereafter 1.2 s pulsa - 
tions were found from Her X-1 (Tananba um et al.|1972 ), 
again with a gradually decreasing period (Ciacconi et al. 
1973). These and other similar sources were soon recog- 


nized as neutron stars accreting from binary companions 

dPringle fc Rees|1972| [Davidson fc Ostriker8l973 

Pethick, . . ^ 


Lamb, 


ethick, & Hines 1973), and the basic theory that was 

erected around them was found to agree approximately 
with their measured spin phenomenology (Rappaport & 
Joss||1977|). The pulses themselves were deduced to be 
caused by channeling of accreting plasma to hotspots 
near the star’s magnetic poles, with the X-rays released 
in the accretion shock subject to a dditional beaming in 
the inferred strong magnetic field (Gnedin & Sunyaev 
1973 Basko & Sunyaevjil975 ). 


'i'wo particularly interesting classes of accreting neu¬ 
tron stars are the millisecond X-ray pulsars and those 
which power relativistic radio jets. Millisecond spin 
periods have been observed for both accretion-powered 
sources, which emit from a magnetically confined accre- 
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tion column, and nuclear-powered sources whose pulses 
are observed during type I X-ray bursts. Their spin fre¬ 
quencies are consistent with a uniform distribution up to 


a cut-off frequency of 730 Hz (Chakrabarty et al. 2003 
Chakrabarty 2008). This is somewhat surprising since 


the cut-ott is well below the critical break-up (mass- 
shedding) frequency estimated for a r ange of internal 
equations of state, about 1230 Hz (e.g. Cook, Shapiro 


& Teukolsky 1994 iLattimer & Hrakash||2004 ). further¬ 
more, for several objects the standard accretion torque 


jjects the standard accretion torque 
theory overpredicts the increase in spin frequency that 
should be measured over t he course of an outburst given 
the observed X-ray flux (Patruno & Watts| 2012), im¬ 
plying either that the accretion torques are less enicient 
than expected or that any torque acting to spin down 
the star is anomalously strong. 

Radio jets have been detected from several binaries 
hosting neutron st ars, most dramatically the non-pulsing 
sources Cir X-1 


d Stew art et al.| |1993|) and Sco X-1 
1981D. Both of these outflows are truly 


(Geldzahler et al.|l98ip. Both of these outflows are truly 


relativistic: energy transfer from the core to the radio 
lobes has been estimated to occur with Lorentz factor 


r > 3 in Sco X-1 (Fomalont, Geldzahler, fc Brad shaw 
2001a) and T > 15 in Cir X-1 ( h'ender et al.||2004[ ). Ex¬ 


tended radio emission has also been associated with sev¬ 
eral accreting millisecond pulsars, which raises the possi¬ 
bility that their jets may be powered by stellar rotation, 
analogous to models of spin-powered black hole jets. Un¬ 
like in black hole systems, there is no evidence for sup- 
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pression of the jet in the soft X-ray state for most ac¬ 
creting neutron stars ( Migliari fc Fender||2006 ), perhaps 
because one or more or the components of the launching 
process is fundamentally different. 

A variety of theoretical models have been proposed 
to describe the interaction between the rotating neutron 
star’s roughly dipolar magnetic field and the highly con¬ 
ducting accretion disk, in tandem with developments in 
the closely relate d pr oblem of p rotostellar accretion; see 
Uzdensky (2004) and Lai (2014) for reviews. Two princi- 
pal states exist, distinguished by whether the disk’s inner 
edge is inside or beyond the corotation radius. 


- 


( 1 ) 


at which the Keplerian and stellar-rotation frequencies 
are equal; M and n are the star’s mass and spin angular 
frequency. If the inner disk is inside r^o it rotates faster 
than the star and accretion along magnetic field lines can 
occur. If the disk is terminated outside corotation a cen¬ 
trifugal barrier prevents accretion, and the rotating field 


the star’s rotational energ-v 

(IDavidson & Ostriker 

1973 

Illarionov & Sunyaev||1975 

Davies & Fringle||1981 

i; this 


is known as the propeller regime. 

The stellar field may penetrate the disk Ghosh , 
Pethick, & Lamb||1977 |Ghosh & Lamb||1978l |l979|), 
resulting in spin-up (spm-downj torques acting on the 
star from field lines entering the disk inside (outside) 
the corotation radius; an accreting star can either gain 
or lose angular momentum, depending on which of the 
spin-up/spin-down torque contributions is larger. In this 
picture the field lines slip through the disk due to the ef¬ 
fective magnetic diffusivity provided by turbulence and 
remain essentially dipolar out to large radii. However, 
even with turbulent diffusivity a thin ther mal disk re¬ 
mains a good conductor, leading Aly (1980) to construct 
a magnetospheric solution in which surface currents suc- 
cessfully exclude the stel lar field from the disk (see also 
Istomin & Haensel 2014). 


Magnetic held lines which couple the star and disk 
are twisted due to the mismatched angular velocities of 
their footpoints; if the diffusivity in the disk is too low 
to allow sufficiently rapid slippage, the resulting build 
up of toroidal field causes the field lines to balloon out¬ 
wards, eventually forming a nearly radial configuration 
in which oppositely directed field lines are separated by 
a current sheet, which may be unstable to reconnection 
(AH 19^ Aly fc Kuijpers I990[ L^mden-Bell fc Body 


994| Uzdensky|2UU2||Ljzderisky, Kmigl, & Litwi^2l)L)2fT 


lOvelace, Romanova, &: Hisnovatyi-Kogan| (|l99o| pro¬ 

posed a model in which the held lines rernam open to 
infinity without reconnecting, allowing magnetohydrody¬ 
namic (MHD) winds to extract angular momentum both 
from t he star and the d isk. In the X-wind configura¬ 


tion of Shu et al. ((r994|), the magnetic field is pinched 
inward and stellar field lines only intersect the disk in a 
narrow annulus around the corotation radius, which ap¬ 
proximately coincides with the disk’s truncation radius. 

The predicted twisting, current sheet formation, and 
reconnection cycle was recovered by the first generation 
of time-depende nt MHD simulations of accretion onto 
magnetic stars (Hayashi, Shibata, & Matsumoto 1996 


Goodson, Winglee, & Bohm 

1997 

Miller & Stone|1997 

'I'he reconnection events anc 

1 subsequently ejected plas 


moids seen in the simulations were associated with the 
X-ray flares and optical jets observed from cla ssical T 


Tauri stars. Kato, Hayashi, & Matsumoto (2004) demon¬ 
strated the acceleration of a transrelativistic jet by the 
twisted magnetic tower crea ted by the neutron star-disk 
system (Lynden-Bell [1996 ). Simulations o f accreting 


stars found solutions with steady funnel flows (Romanova 


et al. 2002 Bessolaz et al. 2008), while those of rapidly 


rotating stars showed the expected propeller-like behav- 


ior (Romanova et al.||2004 Ustyugova et al.|[2006). 


The equations of non-relativistic MHD, used in the 
simulations described above, are satisfactory when 
studying protostars and even neutron stars with peri¬ 
ods of around a second, for which the most important 
interactions between the accreting matter and the stellar 
field occur well within the light cylinder. 


7’LC = TT « 5 --r*, 

il 1 ms 


( 2 ) 


where P and r* are the stellar spin period and radius. 
(In all numerical estimates we use M = 1.4 M© and 
r* = 10 km.) For millisecond pulsars the light cylinder 
is clearly very close to the inner magnetosphere where 
the magnetic field configuration determines the various 
torques the star experiences, and so in principle a rela¬ 
tivistic treatment is required. 

Here we present idealized calculations performed in 
the framework commonly used to describe the extraction 
of rotational energy from an isolated pulsar by an elec- 
tromagnetically dominated plasma wind. We will show 
that, rather than providing a small correction to exist¬ 
ing theory, the increase in the strength of the relativis¬ 
tic pulsar wind due to the star-disk interaction may in 
some cases make this the primary mechanism by which 
spin-down torques are applied to the star. We will dis¬ 
cuss applications to the accretion torques inferred during 
millisecond pulsar outbursts, the spin distribution of the 
millisecond pulsar population, and the production of rel¬ 
ativistic jets in neutron star X-ray binaries. 

We outline in Section the general principles of ac¬ 
cretion onto stars whose internal dipole field is strong 
enough to terminate the Keplerian accretion flow above 
the stellar surface, and describe several models which 
have been advanced to calculate the torque exerted on 
the star by magnetic field lines which intersect the disk. 
In Section we investigate the consequences of the open¬ 
ing of a large fraction of these star-disk coupling field 
lines, and develop a highly simplified model which per¬ 
mits concise algebraic relationships between many of the 
quantities of interest. Sectionj^is devoted to spin equilib¬ 
rium, and the conditions under which it can be achieved 
in our toy model. In Section we compare the behavior 
which would be expected uiider our model to observa¬ 
tions of the spin frequencies, frequency derivatives, and 
jet phenomenology of neutron star X-ray binaries. Fi¬ 
nally we summarize and discuss our results in Section 

2. STANDARD THEORY OF ACCRETION ONTO 
MAGNETIC STARS 

2.1. Basic principles 

A star with a sufficiently strong magnetic field will pre¬ 
vent the direct accretion of matter onto its surface—at 
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the magnetospheric radius the disk is truncated and 
the accreting plasma is directed onto the stellar magnetic 
field lines, along which it funnels down to a ring around 
the magnetic poles, or ring segment if the dipole moment 
is not perpendicular to the disk plane. An estimate of the 
truncation point is given by the Alfven radius where 
the magnetic and kinetic energy densities are equal; it 
can be approximated by equating the magnetic pressure 
from the stellar dipole field to th e ram pressure of gas 
in spherical free-fall from infinity (Davidson & Ostriker 
197^ [Eisner k Lamb||1977D , 


form 


TA = 


( -, (3) 


where ^ is the stellar magnetic dipole moment and M is 
the mass accretion rate. It is standard practice to posit 
the relationship 

Tm = (4) 

in which additional effects, such as from the non- 
spherical disk geometry, are subsumed into the parame- 
ter i < 1. Simulations have generally found ^ ^ 0.5 (e.g. 


Long et al.||2005[ [Bessolaz et al.pOOS Zanni fc Ferreira] 
2013|). Accretion can take place only if < rco- 
Matter accreting from the magnetospheric radius ap¬ 
plies a torque to the star, 


N,co = Ml{r^)^M^yGMr, 


(5) 


where l{r) is the specific angular momentum at radius r, 
given by approximately Keplerian rotation; this torque is 
applied by magnetic field lines twisted by the inflowing 
plasma. Outside the accretion funnel, the stellar mag¬ 
netic field is twisted forward or backward due to drag¬ 
ging by the disk, and any field lines which are open (i.e., 
do not return to the star or intersect the disk) are swept 
back due to the effective i nertia of the electromagnetic 
field (e.g. Bogovalov 1997); together, the disk-coupling 
flux and open fiux apply a torque 


A^field — A^up A^down 


( 6 ) 


to the star, where N^p (A^down) represent the torque con¬ 
tributions which would cause the spin frequency to in¬ 
crease (decrease). The total torque is the sum 

A^tot = A^acc + A^field- (7) 

2.2. Stellar torque from the disk-magnetosphere 
interaction 

Several models have been proposed for estimating the 
torque on accreting pulsars. These are often expressed 
in the form 

A^tot — ‘^(^fast) A^acci (8) 

where n(a;fast) is a dimensionless function and 


n 


^fast — 


f^K(?'m) 


^CO 


3/2 


(9) 


is kn own as the fastness parameter (Eisner & Lamb 
1977); models desc ribed in this way are only applica¬ 
ble when Tm < Tcq. Ghosh & Lamb (1979) proposed the 


^(^fast) — 1.39 


1 ^fast 


4.03 (l-cdfast)®-^^''-0.878 


1 ^fast 

(10) 

although it w as later show n that this model is internally 
inconsistent (Wang 1987). If the magnetic stress com¬ 
municated by the magnetosphere is limited by i ts sus¬ 
ceptib ility to field line opening and reconnection, jWangj 
(19951 argued that the torque should be given by 


^(^fast) — 


7/6-(4/3)a;faat + (l/9)d;Lt ^ 

1 ^fast 


( 11 ) 


Accretion onto millisecond pulsars was addressed by 


(2004), who argued that 


Rappaport, Fregeau, & Spruit 

the disk's inner edge would adjust to be near the coro¬ 
tation radius when the usual estimates, such as equa¬ 
tion (1^, imply that r„j > rco (he., Wfast > !)• In their 
modeb the effective accretion torque is set by the specific 
angular momentum at rco, and the total torque applied 
is given by 


Abot — M \/ GMrco + Afield, 


( 12 ) 


where the contribution from field lines that couple the 
star and disk is 


A^field — 


3r3 


2 2 1 

q -1-2^ 

6 Wfast _ 


9r3 ’ 


a^fast ^ 1 


a^fast ^ 1- 


(13) 


Both this model and equation (11) are constructed us¬ 
ing a magnetic-field pitch distribution in the disk of 
B^/Bz = (11 - Dk) /max(ll, Dk). 

3. OPENING OF PULSAR MAGNETIC FLUX 
BY AN ACCRETION DISK 

3.1. Spin-down of an isolated pulsar 

In the absence of accreting plasma, an isolated star 
rotating about its magnetic axis loses energy at a rate 


,114 


= -Non = 


(14) 


where 'b’open.o is the total m agnetic flux open through 


the light cylinder to infinity ( Gol dreich & Julian 1969_ 
Gruzino^|2005 Gontopoulos||2(Jl)5 ) and Ap is the torque 


applied to the star by the magnetosphere. We will gener¬ 
ally restrict our attention to axisymmetric models, which 
should be a good approximation for stars whose spin and 
magnetic axes are nearly aligned. We take torques which 
increase the stellar rotation rate to be positive, so spin- 
down torques such as A"o are negative. For most accreting 
stars the torque from this electromagnetic wind is small 
in comparison to that due to interaction of the stellar 
magnetic field with the disk, and it is usually ignored. 
However for neutron stars with either large magnetic mo¬ 
ments or high spin frequencies this torque ceases to be 
negligible, and we will argue in particular that Nq sets 
the spin-down torque scale for accreting millisecond pul¬ 
sars. 
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Figure 1. Schematic magnetic field and accretion flow geome¬ 
try for the proposed model, showing the disk outside rm and the 
accretion funnel inside rm (gold region). Most stellar flux which 
would otherwise couple to the disk is forced to open (red lines), 
although a small amount of flux may be trapped near the corota¬ 
tion radius (purple line). Field lines which would be open in the 
magnetosphere of an isolated pulsar of the same spin frequency H 
are indicated in blue, (a) Stellar field lines have opened, forming 
a current layer (cyan dashed line); red arrows indicate the field’s 
poloidal direction on either side of the layer, (b) Relaxed state, in 
which those parts of the open stellar field lines which are rooted in 
the disk have been pushed beyond the light cylinder. 


As a n example, take the simulations of|Ustyugova et ah 


(2006), who measured the spin-down due to the pro¬ 
peller effect on accreting plasma held outside the corota¬ 
tion radius. Scaling their results to a neutron star with 
fi = lO^^Gcm^ rotating at = 5 x 10^ s“^, they report 
a spin-down torque of —3.9 x 10^^ gcm^ s“^. This is ac¬ 
tually slightly weaker than the isolated-pulsar spin-down 
torque for these parameters, Nq = —4.6 x 10^^ g cm^ s“^. 


3.2. Proposed torque & jet power model 

Magnetic field lines which connect the star to the disk 
are forced to open, unless the disk permits azimuthal 
field-line diffusion sufficiently rapid to prevent the ac¬ 
cumulation of toroidal field. Here we focus on the case 
where azimuthal diffusion is not fast enough to inhibit 
opening, and construct a toy model to investigate the 
potential effect the additional open flux may have on 


various properties of the system. 

In the open state, field lines of opposite direction are 
separated by a thin current sheet which will likely be 
unstable to reconnection. The average state of the sys¬ 
tem will then be determined by a competition between 
the rates of differential rotation (i.e., field line twisting) 
and reconnection. When the twisting rate is slow, re¬ 
connection of all the inflated field lines can occur be¬ 
fore the newly reconnected flux has time to expand out¬ 


ward again, leading to cyclical behavior (e.g. Parfrey, Be¬ 


loborodov, & Hui 2012 2013). In the case of accreting 


millisecond pulsars the ditterential rotation rate between 
the star and most points in the disk amounts to hundreds 
of Hz; the conse quent rapid twisting can delay the onset 
of reconnection (jParfrey et al. 2013) and may lead to a 
long-term quasi-steady open state. 

An additional consideration is that, as field lines be¬ 
come increasingly twisted they are driven toward a h- 
nal configuration which is approximately radial in the 
poloidal plane. If the field lines are nearly frozen into a 
conducting equatorial disk they bulge outward at an a n- 
gle of ~ 60° to the vertical (Lynden-Bell & Boily|1994)— 
the field lines are therefore pinched radially inward where 
they enter the disk, and they will slip outward , relative 
to the disk material, on a resistive timescale (|van Bal- 


legooiienl 1989[ Bardou fc HeyvaertsI 1996| [Agapitou &: 
Papaloizou ||2000 [Uzdensky et al.||2002p . If the slippage 
speed is smaller than the radial accretion velocity, the 
segments of the opening field lines connected to the disk 
will move inward; if the slippage speed is greater, the 
field lines will diffuse outward and eventually beyond the 
light cylinder. 

The accretion flow and magnetic field configuration are 
shown in Fig. In Fig. [^a) a strong current layer is 
forming between the two sides of an outwardly inflating 
field line, and the opening field lines are bent forwards (or 
pinched) where they meet the disk; in Fig. Bb) a relaxed 
state has been reached, as the open field lines which are 
now attached solely to the disk have been pushed through 
the light cylinder; the state in (a) may lead to that in 
(b) if the effective magnetic diffusivity in the disk is high 
enough. We have implicitly assumed that the disk’s own 
magnetic field is present only on small scales, where it 
contributes to the turbulent magnetic diffusivity by al¬ 
lowing the stellar field to move through the disk via re¬ 
connection; the disk field may also mediate the coupling 
between the stellar field and the disk. 

Field lines which enter the disk very near the corotation 
radius are only subjected to slow twisting, and may, in 
isolation, be capable of slipping azimuthally fast enough 
to prevent the accumulation of toroidal field. However, 
the magnetic pressure will quickly increase on field lines 
entering well inside r^o, as these field lines experience 
rapid twisting; within a few radians of differential rota¬ 
tion the pressure will be large enough to push outward, 
and open, that small amount of m agnetic flux couplin g 
near the corotation radius (see e.g. Parfrey et al.||2013 ). 
For this reason, the held line crossiiig the equator at the 
corotation radius is, in general, not the same held line as 
would have done so in the disk’s absence. In this way the 
disk-coupling held can be opened (and, depending on the 
accretion velocity, expelled through the light cylinder) all 
the way down to the magnetospheric radius. 

There will be some hux trapped near corotation in 
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the steady state, as illustrated by the magenta field line 
in Fig. [^b); if one considers a sequence of such states 
with decreasing disk magnetic diffusivity, the amount of 
flux trapped in the disk will become infinitesimal as the 
disk approaches a perfect conductor. In this figure, the 
trapped field line entering the disk near corotation will 
have come from near rm if the disk were highly conduct¬ 
ing, and much of the newly opened (red) field lines would 
have initially closed inside Tcq. 

We assume that the open field lines thread a tenuous 
but highly conducting plasma, which for the purpose of 
constructing a simple model we can take to be inertia¬ 
less and dissipationless (i.e., the ideal force-free approx¬ 
imation). Charged particles on the open field lines will 
move radially outward at the Ex B drift velocity, so the 
plasma must be continuously supplied. Radio pulsars 
self-generate this plasma via magnetospheric pair pro¬ 
duction, which may also operate around accreting stars 
if the dense accreting plasma is not able to cross field lines 
onto the open flux, allowing the formation of electrostatic 
vacuum gaps and pair discharge. For a star of a given 
magnetic moment, there is a minimal spin frequency re¬ 
quired for pair production; our sources of interest must 
satisfy this condition, as radio millisecond pulsars are ob¬ 
served with similar magnetic field strengths and rotation 
rates. In an alternative scenario, some of the accreting 
plasma may diffuse onto the open field lines, supplying 
the required charge carriers without the need for pair 
creation. 

The torque on, and power extracted from, the star are 
independent of the source of the wind plasma, and are 
insensitive to its density as long as the open flux region 
is magnetically dominated. The torque and power are 
larger if the open-zone plasma inertia is not negligible, 
because mass-loaded field lines are swept back or twisted 
at a larger angle, increasing the stress applied to the 
stellar surface. The stress is minimized when the only 
inertia coupled to the open field lines is the effective in¬ 
ertia of the electromagnetic field itself. In this sense the 
estimates which follow are conservative, as mass loading 
on the field lines will generally increase the magnitude of 
the spin-down torque or jet power (and reduce the jet’s 
maximum Lorentz factor). On the other hand, the sys¬ 
tem may behave very differently if the accretion rate is 
so large that the disk reaches the star and its entire sur¬ 
face is covered with dense plasma; our estimates are not 
expected to be reliable in this case. 

The amount of stellar magnetic flux that is opened 
following coupling to the disk is difficult to calculate re¬ 
liably, given uncertainties in the general behavior of MRI 
turbulence and the specific complications introduced by 
the stellar field-disk interaction. In the interest of sim¬ 
plicity we therefore subsume the physics that sets the 
amount of open flux into a model parameter, and in¬ 
vestigate the consequences of its taking different values. 
The precise value of C can be expected to depend on the 
conditions in both the accretion disk and the low-density, 
highly magnetized corona. 

If the disk is a good conductor, we can make the esti¬ 
mate that a fraction C ^ 1 of th® stellar field lines that 
intersect the disk are opened. If the field is initially dipo¬ 


lar, the poloidal flux function is 


pd p27T 

tjj{r, d)= / Br r^s'mO dO d(j) = 27r/i 

JO Jo 


sin 9 


(15) 


and at the equator tjj oc 1/r. We therefore estimate that 
when Tni < tlc th® open flux is given by 


/ — I 

V^open — S Y^open,0 5 
'f 'm 


(16) 


where ijjopen,o is the open flux of a n iso lated star at the 
same spin frequency (see equation 14); the conducting 
disk can only increase the open flux and so CcLc/cm > 1- 
The C parameter describes how efficiently the stellar 
magnetic field is opened by the star-disk differential ro¬ 
tation, from C = Tm/cLC (no disk-induced opening) to 
^ = 1 (maximum efficiency or most optimistic case; all 
stellar field lines which intersect the disk are opened). 
The corresponding spin-down torque is 


iVdown.open = CM ^0- 


(17) 


^ig. Qa) results i 


m a 


Finding the torque in this way, by replacing ipopenfi 
ipopen in equation (141, should be a good approximation 
for the fully relaxed^ate of Fig. f^b); th^e poloidal com¬ 
pression of the open flux seen in big. 
slightly larger torque. 

If all of the stellar field is either closed inside rm or 
opened by the disk, the enhanced pulsar wind is the only 
source of spin-down torque, and so ACjown = fVdown,open- 
Likewise, if there are no field lines coupling the star to 
the disk inside corotation, the spin-up torque outside the 
accretion column vanishes, A^up = 0. If the amount of 
flux coupling the star and disk is small (i.e., C ~ 1) we 


can expect A^up Abtcc and A(down ~ ACiowjj 
app roximation we can combine equations 


npcnj, 


in this 
and 


(17) to derive a simple minimal model for the torque, 


Mx/GMtM- r^<r, 


Ntnt — ' 


2 

^ 9 ’ 

rt c 






<rm < r-LC 


rm > rue- 


This model is only appropriate for rapidly rotating stars, 
with spin frequencies greater than roughly a hundred Hz. 
For slower pulsars, such as those with periods of around 
a second, the pulsar wind is much weaker and the spin- 
down torque will be dominated by those magnetic field 
lines linking the star to the disk, even if most of the flux 
is open. 

The pulsar wind may manifest itself as a radio jet 
if it is collimated by external ma gnetic fields or heavy 
plasma, such as a disk wind (e.g. Sulkanen & Lovelace 


1990 Meier, Koide, & Uchida 20(51). As field opening 


adds more flux to the jet, the jet power Lj will be en¬ 
hanced by the same factor as the spin-down torque, 


Lj = CM — 

rm 


Ln. 


(19) 
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Finally, by changing the relationship between spin fre¬ 
quency and spin-down power the disk-induced open¬ 
ing of magnetic flux will necessitate a modified esti¬ 
mate of the pulsar’s magnetic moment. The stan¬ 
dard force-free pulsar lumino sity for general o bliquity 
is Lq = (1 + sin^0)/r^n"‘/c^ (Spitkovsky 20061, where 
9 is the inclination angle between the magnetic and ro¬ 
tational axes, allowing the magnetic moment to be in¬ 
ferred as = /z>c^/[47r^(l -|- sin^d)i^^]. If the disk 
is inside the light cylinder and accretion doesn’t occur 
(rco < ?’m < ?’Lc)i a corrected form should be used; equa¬ 


tion (191 implies 


_ ' m 

— S Mstd 

rLC 

V\ / . x -1/3 

(2GMM^ 


I VC 


-I 7/6 


(1 -I- sixr9)v 


( 20 ) 


Failing to take account of any additional field lines 
opened by interaction with the disk leads to an over¬ 
estimate of the stellar magnetic moment. 


3.3. Numerical estimates 

Using the preceding equations we can estimate the 
spin-down torque for parameters characteristic of accret¬ 
ing millisecond pulsars in outburst. 



( 21 ) 


For the indicated parameters, C = 1; and ^ = 0.5, this 
gives A^down/-^o = 29.3 (reducing v to 250 Hz gives 
Ndov^n/No = 117). 

The jet power can be estimated, using Lj = 
27r7'A^downj as 



which appears to be sufficient to power the most powerful 
neutron star jets (see Section]^; the Eddington accretion 


rate per solar mass of the central star is MEdd,© = 2.3 x 
10“® Mq yr“^, assuming unity radiative efficiency. We 
also find Lj oc r~^ oc we will discuss this in the 

light of observations of the radio-X-ray correlation in jet¬ 


launching neutron star systems in Section 5.3 


3.4. Comparison to other torque models 

In Fig. [^we show the frequency derivative of a star, 
rotating at 500 Hz with a 10^® G cm^ dipole moment, as a 
function of accretion rate, for five cases in which the total 
stellar torme is given respectively by iVacc = My/GMr^ 
(equation^, the three earlier models which include the 


disk-magnetosphere interaction (equations and 

13), and our enhanced pulsar-wind model (equation [T^. 
Here and in Section we relate the torque to a frequency 
derivative using i/ = N/2 ttI and a moment of inertia 


I = 10^® gcm^. 

Discontinuities are present in the v curves predicted 
by our model as we have not attempted to smooth the 
onset of accretion as the magnetospheric radius moves 
inside corotation—in reality these transitions would be 
smeared over some range in accretion rate. The edges in 
Fig. [Kb) serve to indicate the accretion rates at which 


r-m = fco for the three values of ^ shown; for M values 
to the left of these edges the system can be said to be in 
the propeller regime, even though the torque here is due 
to spin-down on open field lines rather than the ejection 
of matter from the disk. 

At high accretion rates, the spin-up torques pred icted 


by previous models—equations 
larger than iVacc due to the inc. 


, and (13)—are 
positive torque 


from stellar field lines entering the disk inside Tco ; m con¬ 
trast, the V predicted by our model is always less than 
the standard estimates, and sometimes significantly so. 
For low accretion rates, where rm <JicQ, our model gives 
stronger spin-down than equation (13). In our model, 
not only can the net torque be negative when the star 
is accreting, but the spin-down can become stronger as 
M increases, as is demonstrated by the ^ = 0.5 and 0.6 
curves in Fig. j^b). The spin-down power on open field 
lines can increase the critical accretion rate at which the 
torque goes to zero, especially for smaller values of 
this will be the subject of the next section. For rapidly 
rotating stars, the open-flux spin-down torque can po¬ 
tentially have a large effect on the star’s spin behavior at 
all but the highest accretion rates. Finally we note that 
the new theory is more sensitive than previous models to 
variations of 


4. SPIN EQUILIBRIUM 

4.1. Balancing accretion and open-flux torques 

Spin equilibrium is obtained when the total torque on 
the star vanishes. In the simple model described above 
(equation[TK) the only spin-up torque is due to the accret¬ 
ing matter and the only spin-down contribution comes 
from the open field lines, and so equilibrium is found 
when A^acc — Adown,open 5 


My/GMr^ 



(23) 


Using equations (|K and (21) we estimate the equilibrium 
spin frequency Vgqim to be 


l/eqlm- 956C ( 1026 G cm3 ) 

( M \ 
\1AMqJ 1^10-10 Moyr-ij 


(24) 


which is only weakly dependent on the mass accretion 
rate. 
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(b) This work, Eq. (18) 
- 5 = 0.5 

- C = 0.6 

- 5 = 0.7 



i''" ■ 


10 -® 
M(MQyr-i) 


Figure 2. Frequency derivative predicted by various models of the total torque Altot, for a star with a 10^® G cm® dipole moment rotating 
at 500 Hz, as a function of mas s ac cretion rate, (a) Earlier models, all with ^ = 0.5, including Ntot = A^acc for rm < rco. (b) Enhanced 
open-flux model, our equation l |18[ l with ^ = 1; the spin-down rate due to the standard isolated-pulsar torque No is shown in the lower 
left. Dashed lines indicate a net spin-down torque. The accretion rates at which rm = rco, and the accretion torque iVacc is switched on, 
correspond to the sharp edges in (b). 


4.2. Maximum steady-state spin frequency for 
accretion-powered pulsars 


Equation (23) can also be rearranged using equa¬ 
tions (pj® ana (fT4| to provide a condition on (or, 
equivalenuy, on 4'j which must be satisfied when spin 
equilibrium is maintained in this way, 


Jjn_ 0 - 1/2 

r-LC ■ 


(25) 


For X-ray pulses to be detectable during outburst the 
accretion flow must be held above the stellar surface by 
the magnetosphere, r^ > f*, implying a maximum self- 
consistent equilibrium spin frequency 




-1 


Hz. 


(26) 


Equilibrium via the above mechanism may still be pos¬ 
sible if the disk reaches the stellar surface, in which case 
this constraint is inapplicable. It is likely however that 
if the accretion rate is significantly larger than the mini¬ 
mal value at which this occurs the open field lines will be 
flooded with heavy plasma at their bases, and the torque 
will differ su bsta ntially from our estimates. 

Equation (26) implies that accretion-powered pulsars 


in instantaneous spin equilibrium are limited to a max¬ 
imum attainable spin frequency, which is independent 
of both the stellar magnetic moment and the accretion 
rate. If C and f only vary within narrow ranges, this 
would enforce a universal cut-off to the spin distribution 
of accretion-powered pulsars, if these objects can be ap¬ 
proximated as being in a steady state at their greatest 
accretion rate. 

Consider a single star with a particular magnetic mo¬ 
ment. There is an accretion rate which places just 



Figure 3. Maximum equilibrium spin fre que ncy versus ^ pa¬ 
rameter, for three values of ^ (see equation |26[ ). Also shown is 
the estimated 730 Hz cut-off to the combined accretion-powered 
nuclear-powered pulsar distribution, and the spin frequency of IGR 
J00291-I-5934, the fastest accretion-powered pulsar (dashed lines). 
The model can match the observational cut-off for reasonable ^ and 
(^, and would set a maximum steady-state frequency independent 
of fi and M. 


above the stellar surface—this is the maximum accre¬ 
tion rate consistent with the star being observable as 
an accretion-powered pulsar, i.e. with pulses formed by 
magnetically confined accretion columns. This accre- 





























tion rate also sets the st ar’s maximum steady-state spin 
frequency via equation (23); if the spin frequency were 
larg er, t he star would spin-down towards i^max (see Sec¬ 
tion 4.4). If M were to increase, the star’s steady-state 
spin frequency would increase, but it would no longe r be 
observed as an accretion-powered pulsar. Equation (26) 
states that the maximum spin frequency found in this 
way is universal for stars of all magnetic moments. 

In Fig. t'max is shown as a function of the model 
parameters. The maximum frequency can be made to 
coincide with the observationally estimated cut-off of ^ 
730 Hz for reasonable values of provided that ( is large 
enough (i.e., that field lines are opened efficiently by the 
disk). 


4.3. Other conditions satisfied by equilibrium solutions 

In this section we describe several additional relation¬ 
ships which must hold for a spin equilibrium solution to 
be physically valid, and show that these conditions are 
satisfied by the observed systems given reasonable ranges 
of the model parameters. 

The star must be accreting for spin equiilibrium to 
be possible, and so we can demand that the magneto- 
spheric radius is inside the corotation radius, r^ < r^o', 
equations 0 and ( |25[ ) then imply that in equilibrium 
the star must be spinning faster than some minimum 
frequency, i^eqim > t'min, where 


= 8160 


M 

AMq 


Hz. 


(27) 


The condition that r^ < rco in spin equilibrium can 
also be understood as requiring that r^ lies inside rm^max, 
the largest possible corotation radius consistent with 
equilibrium at a given ^ and which is 

^m.max ~ ^co(^min) ~ C ^ (^^) 

where rg = 2GM/c^ is the star’s Schwarzschild radius. 
Demanding r^ < rm,max in turn implies that M > Mmm, 


= 1.81 X 

Although the very strong dependence on the C and ^ 
parameters makes reliance on this quantity dubious, we 
can see that for, say, C = 0-9 and ^ = 0.7 one finds 
.^min = 3.6 X 10“^^ Mq yr“^, distinctly lower than the 
observed outburst accretion rates for the millisecond pul¬ 
sars; as ^ is reduced further, Mmm becomes small very 
rapidly. 

Clearly, we require r-min < t'max; since t-min/i^max « 
2.42 for our standard stellar mass and radius, this 
provides a maximum self-consistent value of ^ for a given 
^ in spin equilibrium, 

^ < 0.881 (30) 

If we wish to explain the observed accretion-powered 
millisecond-pulsar spin distribution as consisting of sys¬ 
tems in approximate equilibrium, we can demand i^max > 



Figure 4. Equilibrium spin period versus magnetospheric radius 
for the standard model (black curve) and our enhanced pulsar wind 
model (colored lines). The equilibrium described here can only hold 
for rm < rm,max, where rm,max is given by the intersection of a 
particular Feqim(C?0 with the rm = rco curve. 


730 Hz, implying a minimum 


e > 0.645 


V730 Hz/ 


(31) 


Note that for zzmax = 730 Hz the ratio of minimum to 
maximum allowable f is comfortably smaller than unity, 
independent of C. 

Smaller values of ^ are required to explain slower equi¬ 
librium frequencies in this framework, because for slower 
spin one must open more field lines for a, given accretion 
rate (i.e. N^ovjn oc see equation 21). It is not clear 

at what frequency the enhanced pulsar wind would cease 
to be the dominant source of spin-down torque, but using 
equation (27) we can estimate an alternative maximum 

The very weak sensitivity to lymin implies that ^ ~ 0.5 
would provide a valid equilibrium solution even for slow 
rotators. 


Equations ( 30 - 32 ) imply that ^ varying within the rea¬ 
sonable rangeV.o ^ ^ < 0.8 would produce valid equi¬ 
librium solutions for the entire millisecond pulsar distri¬ 
bution. Values of ^ covering much of this range (in par¬ 
ticular 0.55 < f < 0.72) were measured in a numerical 
study of the depend ence of on /i and M by Kulkarni 
& Romanoval (|2013|). 


4.4. Relating equilibrium rotation rate to 
magnetospheric radius 


We could alternatively think of equation ( [^ as ex¬ 
pressing the equilibrium period Peqim predicted by the 
model. 


2^/^7r 2 ,,_ 7/2 

/^qlm = - 


(33) 


c 
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Figure 5. Fastness parameter cjfast = f^/^K(^m) versus spin 
frequency, at spin equilibrium. For a valid equilibrium solution 
^fast,eqim Hiust be less than unity. 


here Peqim oc r„i- In th e standard p icture of spin equi¬ 
librium, Tm « r-co (e.g. 


Wang 19951 and so Peqim.std 


2'k{GM) . These two scenarios are illustrated in 

Fig-i In this diagram a star spinning up at constant 
accretion rate moves vertically downward, assuming ^ 
remains constant. 

If Tm < rni.max (sBC equation 28), the star will eventu¬ 
ally reach spin equilibrium at trie period given by equa¬ 
tion (33). The torque at any stage is determined by the 
competition between iVacc, which is not a function of 
the spin period, and Vdown oc P~^. When P > Peqim, 
Vacc + Vdown > 0 and the star spins up, decreasing P 
and increasing the strength of the spin-down torque un¬ 
til an equilibrium is reached. Likewise when P < Peqim, 
Vacc + Vdown < 0 and the star’s period increases, de¬ 
creasing the spin-down strength until Vacc + Vdown = 0. 
In this way equation (33) describes an attractor for P. 

If rm > r-ni^max, the corotation radius, which moves in¬ 
ward as the star is spun up, reaches before the star is 
rotating fast enough to establish spin equilibrium using 
the open flux set by r^. In this case the star’s rotation 
will hover around Feqim.std: accretion is shut off when 
fco < fm, the star spins down, r^n increases, until accre- 


tion and spin-up restart (see e.g. Spruit & Taam 1993 


D’Angelo & Spruit 2010). This long-term averaged (as 


opposed to instantaneous) spin equilibrium may operate 
in systems hosting neutron stars with longer periods. 

Fig. 1^ suggests that accretion torques may be instan- 
taneou^y balanced by spin-down torques on open field 
lines for some of the most rapidly rotating stars, in the 
most extreme cases preventing spin-up to sub-millisecond 
periodS; _ 

The |Ghosh & Lamb| (|1979|) and |Wang| (|1995| torque 
models, equations ([tOp and ([lT|) , predict fastness param¬ 
eters at spin equilibrium (often referred to as the critical 
fastness parameter) of Wfast.eqim = 0.35 and 0.95 respec¬ 
tively, independent of r'eqim- Using equations ([^ and 


(25) we can express this quantity in our model as 


Wfast.eqlm = 2 "V4^-3^21/4^3/2 (GMU)"'/" . (34) 

In Fig. Wfast.eqim IS plotted against z/gqini for several 
combin^ons of C and A valid spin equilibrium so¬ 
lution requires rm < Vco and hence Wfast,eqim < 1; the 
locations where Wfast.eqim = 1 therefore coincide with the 
minimum allowable spin frequencies of equation (27). As 


the equilibrium frequency increases, the larger spin-down 
power must be balanced by greater mass accretion and 
hence a decreasing magnetospheric radius, which is re¬ 
flected in Wfast.eqim being a decreasing function of t'eqim- 


5. COMPARISON TO OBSERVATIONS 
In this section we express magnetic moments as /i = 
fi 26 X 10^® G cm3, accretion rates as M = M_io x 
MQyT~^, torques as N = N 32 x 103^gcm^s“^, 
and frequency derivatives as z> = z>_ 14 X 10-14 Hzs-b 
Where necessary we convert reported surface magnetic 
field estimates B to dipole moments using p, ~ Br^; when 
spin-down frequency derivatives in quiescence are avail¬ 
able we calculate /r directly using the force-free aligned 
rotator formula. 


5.1. Spin frequency distribution of millisecond 
X-ray pulsars 

The frequency distribution of rapidly rotating pulsars 
has been puzzling at least since White & Zhang (1997) 
reported an apparent correlation between magnetic held 
strength and accretion rate, required to explain the exis¬ 
tence of pulsars of similar, presumably equilibrium, peri- 
ods at w i dely v arying accretion rates. This question led 
Bildsten (1998) to suggest that the strongly frequency- 
dependent onset of spin-down via gravitational radiation 
enforces a maximum cut-off rotation rate; alternatively, 
the spin-up torques may be reduced at high accretion 
rates due to the formation of a radiation-pressure sup - 
ported sub-Keplerian inner disk (Andersson et al.|2005). 
While the discovery of more objects has obviated con¬ 
cerns over a magnetic field-X-ray luminos ity correlation 
(Patruno et al. 2012 Papitto et al. 2014), the problem 
of explaining the lack ol pulsars above ~ 700 Hz remains 
unresolved. 

We described in Section |3] how interaction with a disk 
leads to increased open stellar magnetic flux, allowing 
a state where the spin-up accretion torque is instanta¬ 
neously balanced by angular momentum loss through the 
pulsar wind, as outlined in Section]^ For magnetic mo¬ 
ments and accretion rates consistent with observations, 
/i 26 ~ 1 and M_io ~ 1-50, and theoretically reasonable 
model parameters C 1 and ^ ^ 0.5-0.8, the model can 
explain the spin frequencies of the observed millisecond 
X-ray pulsars as being set by the termination of spin-up 
by the sta r en tering th e spin equilibrium state given by 
equations (23) and (24). 

In particular, the model predicts a maximum steady- 
state spin frequency for sources which appear as 
accretion-powered pulsars, as set by equation ( |26[ ); see 
Fig. This cut-off is independent of the accretion rate 
and the star’s magnetic field, and would therefore be uni¬ 
versal if all magnetosphere-disk interactions can be de¬ 
scribed by similar f and C; it can be reconciled with the 
observationally estimated ~ 730 Hz maximum frequency 
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for the same reasonable range of model parameters as 
above. 

A model in which the primary limitation on an ac¬ 
creting star’s rotation rate is due to magnetically ap¬ 
plied torques would accord with observations which sug¬ 
gest that the nuclear-powered pulsars (burst oscillation 
sources) have a faster spin distribution than both the 
accret i on-po wered and radio millisecond pulsars (Papitto 
et al. 2014). The nuclear-powered pulsars generally 
have a weaker surface magnetic field than the accretion- 
powered pulsars, giving a smaller spin-down torque for 
the same accretion rate at all stages of an outburst, and 
possibly periods of little or no spin-down torque when 
the stellar surface is flooded at the highest accretion 
rates. These objects are not stri ctly subject to the high- 
frequency cut-off of Section |4.2[ but their maximum at¬ 
tainable spin frequency can oe expected to be similar 
to that of the accretion-powered sources, as appears to 
be the case, if the magnetic field strengths of the two 
populations don’t differ too substantially. The frequen¬ 
cies of the nascent radio millisecond pulsars may have 
been systematically reduced by the spin-down mecha¬ 
nism described here during the terminal, low accretion 
rate, p hases of Roch e-lobe overflow, in a manner analo¬ 
gous to Tauris (2012), who used a variant of equation (131 
for the propeller-regime spin-down torque. 

Another possible scenario is that the star’s increasing 
spin frequency may not be able to reach the equilibrium 
frequency commensurate with the outburst’s maximum 
accretion rate before the end of an outburst, and the 
frequency increase sustained by the star may be approx¬ 
imately balanced by strong net spin-down due to the en¬ 
hanced pulsar wind as the accretion rate declines; see 
Fig.i In this way the stellar spin frequency may oscil¬ 
late in a sawtooth pattern over the outburst-recurrence 
timescale while being in a kind of long-term spin equilib¬ 
rium. 

5.2. Spin evolution of individual sources 

For several sources, frequency measurements during 
multiple accretion episodes allow the spin-down rate dur¬ 
ing X-ray quiescence to be found, with which one can 
estimate the star’s magnetic moment using the force-free 
spin-down formula. This in turn permits tests of the 
torque models during, or integrated over, a single accre¬ 
tion outburst. This is complicated by the presence of 
strong timing noise in the observations, making the cal¬ 
culation of reliable v's much more problematic than for 
isolated pulsars, and the nonlinear response of both the 
accretion and spin-down torques to the accretion rate. 
Despite these difficulties, a concern has been voiced that 
some of these objects do not show as strong a spin-up 
during outbursts as the X-ray luminosity would imply 
under the standard models (e.g. Patruno fc Watts|2012 ); 
in some cases no change in frequency was discernable 
despite the inference of a large accretion rate. Here we 
briefly show that the anomalous spin-down torque re¬ 
quired is ap prox imately consistent with the torque model 
of equation (181 for reasonable choices of the various pa¬ 
rameters, and describe similar observations of other sys¬ 
tems which do not require additional spin-down. We take 
Ntot = A^acc = as the standard torque model 

at high accretion rates (see Fig. [^. 


An RXTE study of the pulse frequency during five 
outbursts of the v = 401 Hz pulsar SAX J1808.4-3658 
by Hartman et al. (2008 2009) has revealed a long-term 
spin-down rate of r'_i 4 « —y (jj 55 This spin-down is con¬ 
sistent with being constant in time, and is too large to be 
due to torques during the outbursts given their duration 
and the single-outburst upper limit of |f'-i 4 | < 2.5, and 
so is assigned to spin-down during X-ray quiesence; this 
implies a dipole moment of p ,26 ~ 0.76 (without torque 
enhancement). 

Using an average outb urst accretion rate of M_io « 5 
(Haskell & Patruno 20111, the spin-up torque in the stan¬ 
dard theory would be Aacc ,32 ~ implying a fre¬ 

quency derivative v-u « 9.4v^, well above the observa¬ 
tional upper limit. The spin-down torque, equation ([TtI), 
would be A^down ,32 ~ 1-4(C/C)^ so the combined 
model, equation (18), predicts !>_i 4 « 9.4v^—2.2(((/^)^. 
For the frequency derivative in the model to obey the ob¬ 
servational constraint (i.e., f'_i 4 < 2.5), one must there¬ 
fore have f < [0.65, 0.61, 0.55] for models in which f = 
[1.0, 0.9, (1.8] respectively; note that > r* for these 
values of As these values of f are reasonable, given 
previous theoretical and numerical studies, the enhanced 
open-flux model can explain the undetectable spin-up of 
SAX J1808.4-3658 if C, is sufficiently close to unity. 

The long-term mass transfer rate in the J1808 system is 
expected to be approximately 10~^^ (Bildsten 

& Chakrabarty 2001 Heinke et al.||2009 ); for tnis value 


the corrected rhagnetic moment given by equation (20) 
would be /4corr,26 ~ 0.15, which is inconsistent with the 
detection of pulsat i ons at high accretion r ates (Psaltis & 


Chakrabarty 1999 Hartman et al. 2008). This may be 
resolved by a more detailed model ot the build up of the 
disk between outbursts. 

Observations of the 2003 outburst of the 314 Hz pul¬ 
sar XTE J1814-338 placed an upper limit of If'- 14 ] < 1. 5 
on the frequency derivative (Haskell & Patruno 2011). 
As only one outburst has been detected by an instru- 
ment with sensitive timing capability, no estimate can 
be made of the quiescent spin-down rate and hence mag¬ 
netic moment. Usi ng the average outburst accretion 
rate of M_io ~ 2 (Haskell & Patruno 2011) one de¬ 
rives, in the standard torque model, an expected spin- 

up rate of v-ia, « 4.7/i2g^v^, well above the obser¬ 
vational constraint for /i 26 ~ 1- Including spin-down 
torques on open field lines one finds i>-\A « 4.7 — 

1.3/i2g^(C/'?)^! giving, for ^26 = 1, scenarios satisfying 
the observational constraints (i.e., |f'-i 4 | < 1.5) when 
^ < [0.72,0.67,0.61,0.56] for C = [1.0,0.9,0.8,0.7] re¬ 
spectively; all of these ^ result in > r* as ta = 

28.92 /i 2 g^ km for this accretion rate. Note that this 
source can be in spin equilibrium during such an out¬ 
burst given spin-down on open flux alone: = 0 for 

^ = [0.59,0.54,0.49] using = [1.0, 0.9,0.8]. We conclude 
that an enhanced open-field torque may be responsible 
for the weak or absent spin-up of XTE J1814-338. 

Two systems for which accretion-induced spin-up has 
been measured do not appear to require enhanced brak¬ 
ing torque. The 2002 outburst of the 435 Hz pulsar XTE 
J1751-305 had a peak X-ray luminosity of Lx « 2.7 X 


10^"^ ergs ^ for an estimated distance of 8.5 kpc (Mark- 
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ward! et al.|2002 Gierlinski fc Poutanen|2005 ), implying 
a mass accretion rate of M « 2.3 x 10“® Mq yi~^ « 
-^Edd, 0 - The peak fr equency derivative h as been re¬ 
ported as v -14 « 56 by Papitto et al. (2008). The long¬ 
term spin-down rate of ;>-i 4 » —0.55 indicates a mag¬ 
netic moment of /i 26 ~ 2.1 (Riggio et al. 2011). The 


above values give an outburst peak frequency deriva¬ 
tive, due to the spin-up torque component A^acc alone, 
of f'-id « 47 v^, smaller than the measured value, and 
therefore additional spin-down torque is not called for. 
A large spin-up during o utburst has also b een ob¬ 


served for IGR J00291-I-5934 (Falanga et al. 2005), at 599 
Hz the fastest rotating accretion-power ed pulsar. The 


—0.41 (|Hartman et al. 


2011| Papitto et al. 2011 ) gives ^26 ~ 1-1. |Patruno 
12010) estimated the peak accretion rate of the 2004 out- 
burst to be M_io ~ 20 and measured a spin-up rate of 
z>_i 4 « 51. The accretion spin-up torque alone would 
give iz-id « 35-\/J for the above parameters, or i>_i4 « 42 
if the torque were applied from the corotation radius. No 
additional spin-down torque component is required to ex¬ 
plain the timing behavior of this system in outburst. 

5.3. Jets from neutron star X-ray binaries 

The radio emission of Gir X-1 and Sco X-1 has been 
resolved into clearly defined jets. The brightening of the 
Gir X-1 radio lobes following an X-ray flare implies en¬ 
ergy transfer at Lorentz factor P > 15, and minimal- 
energy arguments for powering the exte nded emission 


requi re a jet power Lj > 10^® ergs ^ (|Fender et al. 


2004). For Sco X-1, energy transfer along the jet occurs 
at i' > 3, and equipartition magnetic fields of ^ 0.3- 
1 G and a radiative lifetime of ^ 30 minutes ag ain de¬ 
mand a jet power in excess of a few 10^® ergs“^ (Foma-, 
lont et aL]|2001a|b ). These objects typically accrete at a 
large traction ot the Eddington rate, and are categorized 
as Z sources in the spectral an d timing classification of 
Hasinger & van der Klis (19891. Both may be rapidly ro- 
tating, as suggested by separations of a few hundre d Hz 
between their twin kHz quasi-periodic osc illations (van 
der Klis et al.|1996 iBoutloukos et al.|2006 ), although see 
Abramowicz et al. (|20U3 ) tor an alternative viewpoint on 
the kHz GPUs, it has been proposed that the Z sources 
have relatively strong magnetic fields; dipole moments 
in excess of 1 0^^ G cm^ have been e stimated for Gyg X-2 
(|Focke||1996|) and XTE J1701-462 (jDing et al.||20ii). 

It the magnetic moments of Gir X-1 and Sco X-l are 
sufficiently large, greater than about 10^® Gcm^, the re¬ 
quired jet powers can be supplied by the rotating open 
stellar field lines whose number has been increased by 
interaction with the surrounding disk. Eor p ,26 = 1) 
v = 300 Hz, and mass accretion at half the Eddington 
rate (i.e. M = O.5[M/M0]MEdd,0) equation (22) gives 

35/ 


Lj = 4.6 X 10^'^(C/O^ ergs“^ C/^ = 1 this is two 
orders of magnitude larger than the spin-down power of 
an isolated pulsar of the same dipole moment and spin). 

When the star’s spin and magnetic axes are nearly 
aligned, the center of the open flux bundle will be mag¬ 
netically shielded from the accreting matter, allowing 
the formation of a Poynting flux-dominated jet with low 
baryon loading and hence high terminal Lorentz fac¬ 
tor; a similar effect was observ ed in the non-relativist ic 
propeller-regime simulations of Romanova et al. (2009). 


Several accreting millisecond pulsars, including SAX 
J1808.4-3658 and IGR J00291-I-5934, are associated with 
radio emission which is attributed to a jet; they have gen¬ 
erally been found to have stronger radio emis sion than 
the atoll sources at the same X-ray luminosity (Migliari 


Miller-Jones, & Russell 2011). This could be explained 
in the open stellar field framework by their having gen¬ 
erally higher spin frequencies or stronger magnetic fields. 
The latter may be the case, as accretion-powered pulsars 
need relatively strong fields to confine the accretion col¬ 
umn and atolls are thought to be weak-field objects; for 
example the magnetic moment of the atoll source Aql X- 
1 has been e s timated to be le ss than 10^® G cm^ (Zhang 


et al. 1998a Campana 2000). Analyses of the relative 
rankings of sources’ rotation rates and radio emission 
have found some evidence for higher spin frequency be¬ 
ing associated with higher radio luminosity, as would be 
expected in a model wher e the jet power ultimately de¬ 
rives from stellar rotation (Migliari et al. 2011 King et al. 


2013). 

The most basic model of jets pow ered by the en¬ 
hanced stellar open flux, equation (221, predicts a com¬ 


paratively weak scaling of jet power with accretion rate 
Lj (X This is similar to what is seen in Aql X-1, 

a rapidly ro tating star with 54 9 Hz type I X-ray burst 
oscillations ( Zhang et al.|1998b ), whose radio luminosity 
Lr scales with X-ray luminosity as Lr oc on the 

other hand , 4U 1728-34 shows a s tronger dependence, 
Lr oc ( Migliari fc Eendei]|2006 ), which is con sistent 
with Lj cx Af in t he scale-invariant disk-jet model (Heinz 
fc Sunyaevj2003 1. In the scenario we describe here, 4/Y is 


a lower limit to the scaling power with the accretion rate 
since we base the jet power on the spin-down luminos¬ 
ity of the standard isolated pulsar wind solution, which 
is the lowest spin-down state for a given open magnetic 
flux. If the disk is thick the poloidal flux will be pushed 
away from the equator, increasing the jet power; as M 
increases, the magnetospheric radius moves inward and 
the amount of poloidal field line displacement increases, 
which corresponds to the jet power scaling somewhat 
more strongly with M. 

Neutron star X-ray binaries have generally shown less 
tendency than those wi th black holes to suppres s the jet 
in the soft X-ray state ( Migliari fc Fend^|2006 ). In the 
model discussed here, jet suppression would be expected 
when the accretion rate increases to the point that the 
disk reaches the stellar surface and all of the field lines 
are inundated with accreting matter. 

Quenching has been obser ved in Aql X-1 at accretion 


rates above about O.lMFUri (Tudose et al. 2009 


Miller- 


Jones et al. 

2010 

). At this accretion rate the Alfven 

radius ta ~ 

lUkm lor a dipole moment ot fi 26 = ^-^5 so 


if the field is relatively weak, fj. 2 e ^ 0-7, flooding may 
be responsible for the jet suppression, depending on the 
effective ^ = r^jr^. Neutron stars with stronger fields 
may be immune to this type of quenching; for > 2 
the Alfven radius is above 20 km for accretion at the Ed¬ 
dington rate, and so the jet should persist even at the 
highest accretion rates unless ^ is unexpectedly small. 
This mechanism would explain why only some binary 
systems, those with stars whose magnetic moments are 
smaller than a critical value of roughly a few 10^® G cm®, 
are subject to jet suppression in the high accretion rate. 
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soft X-ray state. This contrasts with jet models based 
solely on accretion disk physics, similar to those em¬ 
ployed for accreting black holes, where quenching may 
be tied to a transition from a thick to a thin disk at 
roughly O.lMEdd; in that case all systems should show 
jet suppression in the luminous soft state. 

Several points in the above discussion would be clar¬ 
ified by reliable measurements of the dipole moments 
of jet-associated neutron stars, particularly as the pro¬ 
posed difference in field strength between the Z and 
atoll sources has been challenged by observations of XTE 
J1701-462, which has been observed to displa y both Z- 


and a toll-like behavior during a single outburst (Lin et al. 


20091. 


6. DISCUSSION & CONCLUSION 

The rotational energy of rapidly spinning neutron stars 
can be efSciently extracted in a powerful Poynting flux- 
dominated wind, even for stars with comparatively low 
magnetic dipole moments, if magnetic field lines which 
would otherwise close within the light cylinder can be 
forced to open to infinity; the strength of the usual elec¬ 
tromagnetic pulsar wind is effectively boosted by a large 
factor. A natural means of opening these field lines is 
coupling to a surrounding accretion disk—differential ro¬ 
tation between the footpoints of a field line linking the 
stellar surface to a Keplerian disk causes toroidal mag¬ 
netic field to accumulate and the field line to inflate out¬ 
ward, eventually achieving an open-field configuration. 

We have presented a simplified minimal model of the 
torque on the star and the power in the electromagnetic 
wind when a fraction of the total stellar flux which would 
otherwise link to the disk is opened. This wind may 
be collimated by external gas or magnetic pressure, or 
by mass loading on its outermost magnetic field lines, 
creating a jet. We have applied our torque and jet power 
estimates to three problems: 

1. The millisecond X-ray pulsar spin-frequency cut-off. 
We find that at high spin frequency, and moder¬ 
ate accretion rates and dipole moments, the spin- 
up torque due to accretion can be balanced by the 
spin-down torque from the open field lines alone, giv¬ 
ing spin equilibrium and preventing spin-up to sub¬ 
millisecond periods. Unlike the proposed action of 
gravitational radiation, in this case the torque does 
not scale as a high power of the frequency and so would 
not enforce a nearly universal maximum spin. Never¬ 
theless, the model predicts a maximum steady-state 
spin frequency for neutron stars which are observable 
as accretion-powered pulsars, a limit which is inde¬ 
pendent of accretion rate and stellar magnetic field 
strength. For individual sources below the cut-off, the 
equilibrium spin frequency depends on parameters like 
the accretion rate and the stellar mass and magnetic 
field, which appears to be generally consistent with the 
observed flat spin distribution of accreting millisecond 
pulsars. 

2. Weak or absent spin-up of millisecond pulsars in 
outburst. Accretion episodes of two sources, SAX 
J1808.4-3658 and XTE J1814-338, resulted in unde¬ 
tectable frequency increases, despite X-ray emission 
which would lead one, using standard accretion torque 


theory, to expect spin-ups well above the observational 
limits. We have argued that the enhanced open flux 
spin-down torque could have been strong enough to 
prevent measurable spin-up in these systems. 

3. Relativistic jets from neutron star binaries. A model 
in which the jet is powered by stellar rotation and 
mediated by enhanced open stellar flux can account 
for the strengths of the Cir X-1 and Sco X-1 jets, for 
the observed accretion rates and reasonable magnetic 
dipole moments and spin frequencies. The field lines 
at the center of the jet would not connect to the disk 
and would be shielded from contaminating baryons, 
which may result in the observed high jet Lorentz fac¬ 
tors. This model also appears to be consistent with 
the higher radio luminosity of the accreting millisec¬ 
ond pulsars compared to the atoll sources, and may ex¬ 
plain the jet suppression in Aql X-1 at high accretion 
rates and the lack of such quenching in other sources. 


The problem of accreting magnetic stars is clearly 
highly complicated, and the idealized picture described 
here can only possibly tell a small part of the story. Our 
basic conclusion is that the open stellar magnetic flux 
could make an important contribution to the behavior of 
many neutron star systems if its share of the total stel¬ 
lar flux can be increased by the presence of surrounding 
plasma. We will present results of numerical simulations 
addressing various aspects of this problem in the context 
of our model in a future paper (K. Parfrey et al. 2016, 
in preparation). 

For simplicity we have assumed that the stellar mag¬ 
netic field and the accretion disk are in a stationary con¬ 
figuration, determined by the accretion rate which can 
be taken as constant; this is only a rough approximation, 
but may often be acceptable, as the timescale for relax¬ 
ation of the magnetosphere to the steady state will be 
related to the spin period of the star, which is very short 
for the objects of interest. There are two cases where 
time dependence may be particularly relevant. Firstly, 
in the early stages of accretion outbursts, when the inflat¬ 
ing stellar magnetic field is competing with the increas¬ 
ing accretion rate and more flux may couple the star and 
disk. In this case the spin-up torque contribution will 
be higher and the spin-down contribution lower than in 
our estimates, leading to stronger spin-up; this may be 
responsible for the large frequency derivatives observed 
over outbursts of XTE J1751-305 and IGR J00291-f5934. 
Secondly, increases in the accretion rate or changes in the 
disk structure coincident with X-ray binary state transi¬ 
tions may trigger the sudden opening of additional mag¬ 
netic field lines, leading to reconnection and the emis¬ 
sion of large plasmoids before the steady state is reached; 
this may be related to observations of powerful transient 
jets a ssociated with state transitions (Migliari fc Fendei] 


20061 . 

Magnetic flux connecting the star to disk plasma 
with lower angular velocity applies a stronger spin-down 
torque to the star than the same flux opened to infinity 
and contributing to spin-down via the enhanced pulsar 
wind, if the disk’s resistivity is low and field lines do 
not easily slip through it. That open field lines ca n gen¬ 
erate less spin-d own torque than closed ones led |Matt| 


& Pudritz (20051 to worry that a largely open configura- 
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tion would have difficulty preventing accreting protostars 
from spinning up to unobserved angular frequencies. We 
also propose that there will generally be little flux cou¬ 
pling the star and disk, however in the case of rapidly 
rotating neutron stars the combination of high spin and 
relatively strong magnetic fields means the open flux 
can nevertheless dominate the system energetically. For 
these objects external mass loading is unnecessary be¬ 
cause the rotationally induced voltages are large enough 
for magnetospheric pair creation to supply the charge 
carriers required by the electromagnetic wind, as mani¬ 
fested in the coherent radio emission of isolated pulsars 
of similar spins and magnetic moments; additional iner¬ 
tia on the open field lines would however increase a star’s 
spin-down power. 

The picture of neutron star jets proposed here is related 
to Blandford-Znajek-type models for black hole jets, in 
that in both cases the jet power is extracted from the ro¬ 
tational energy of the compact object by helical magnetic 
field lines; it makes little qualitative difference whether 
the power transfer from the compact object to the field 
is mediated by a rotating space-time or a rotating solid 
surface. This outlook contrasts with magnetocentrifugal 
jet models, in which the jet is launched by field lines 
attached to the inner accretion disk; it has been argued 
that to produce a jet in this way the stellar magnetic mo¬ 
ment must be small e nough that the disk extends down 
to the stellar surface ( Mass! fc Kaufman Bernad6||2008 ). 
Both mechanisms may operate in the same system at dif- 
ferent accretion rates or even concurrently—in particular 
a disk-launched jet may persist when the stellar-field jet 
is quenched at high accretion rates. 

The lack o f observed jets fro m strong-held pulsars in 
binaries (e.g. Fender et al. 19971 may also be explained in 
this framework—here jet power scales like Lj oc 
and strong-held pulsars’ high dipole moments are out¬ 
weighed by their slower rotation, to which the jet power is 
more sensitive; the jet power from a star with B ~ 10® G 
and V ^ 500 Hz is about 100 times larger than from one 
with B ~ 10^^ G and v ^ 1 Hz, for the same accretion 
rate and 

An interesting new class of objects is the transitional 
millisecond pulsars: neutron stars which switch between 
X-ray pulsar and ra dio pulsar phases, over timescales as 
short as weeks (e.g. Archibald et al.||2009 Papitto et al. 


2013| [Bogdanov et al.||2014|). (Jt these sources, only tor 
PSR J 1023-1-0038 has a magnetic moment been estimated 


, magnetic i 

by measuring a freque ncy derivative in the radio phase 
( Archibald et al.||2013 ), and there are no constraints on 
frequency derivatives during accretion episodes. The in¬ 
crease of the apparent dipole moment inferred from spin- 
down measurements due to the opening of additional stel¬ 
lar flux (see equation 20) may be particularly important 
for the transitional pulsars, whose pulsed radio emission 
may resume before the remnant disk is cleared from the 
light cylinder. It is hoped that the unique possibilities of 
these systems may prove useful for disentangling the re¬ 
lationships between accretion rate, stellar magnetic field, 
spin-up/spin-down torques, and jet launching. 
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